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A five-qubit codeword stabilized quantum code is implemented in a seven-qubit 
system using nuclear magnetic resonance (NMR) . Our experiment implements a good 
nonadditive quantum code which encodes a larger Hilbert space than any stabilizer 
code with the same length and capable of correcting the same kind of errors. The 
experimentally measured quantum coherence is shown to be robust against artificially 
introduced errors, benchmarking the success in implementing the quantum error 
correction code. Given the typical decoherence time of the system, our experiment 
illustrates the ability of coherent control to implement complex quantum circuits for 
demonstrating interesting results in spin qubits for quantum computing. 

PACS numbers: 03.67.Pp, 03.67.Lx, 03.65.Wj 



INTRODUCTION 



Quantum computers are promising to solve certain problems faster than classical com- 
puters [l|. The power of quantum computing relies on the coherence of quantum states. In 
implementation, however, the quantum devices are subject to errors, from inevitable cou- 
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pling to the uncontrollable environment, or from other mechanisms such as imperfection 
in controlled operations. The errors damage the coherence, and consequently can reduce 
the computational ability of quantum computers. In order to protect quantum coherence, 
schemes of quantum error correction and fault-tolerant quantum computation have been 
developed jl-7|. Those schemes have greatly improved the long-term prospects for quantum 
computation technology. 

A quantum error correcting code (QECC) protects a /T- dimensional Hilbert space (the 
code space) by encoding it into an n-qubit system, and is usually denoted by parameters 
((n, K,d)), where d is called the distance of the code js]. This n-qubit system is used in the 
process of quantum computing and hence subject to errors. A code of distance d is capable 
of correcting d — 1 erasure errors (i.e. loss of qubits at up to (i — 1 known positions) or 
t = [^^J arbitrary errors (i.e. arbitrary errors on t qubits at unknown positions). At the 
end of the computation, the quantum code can be decoded to recover the quantum state in 
the original i^- dimensional Hilbert space. 

In practice, one would always hope for a "good" QECC where more information is pro- 
tected (larger K), while less physical resources are used (smaller n) plus more errors can 
be corrected (larger d). There are trade-offs among these three parameters, and one can 
readily develop upper bounds and lower bounds for the third parameter if two of them are 
fixed [7j . With increasing length n of the codes, in most of the cases there is a gap between 
these upper and lower bounds. Therefore, given two fixed parameters among the three, to 
find a code with the best possible value of the third parameter is one of the most important 
topics in studying the theory of QECC. 

Stabilizer codes, also known as additive codes, form an important class of QECCs, de- 
veloped independently in {g, 7| in the late 1990s. The construction of these codes is based 
on a simple method using Abelian groups, where the code dimension K of stabilizer codes 
is always a power of two, that is, K = 2^ for some integer k. Stabilizer codes include the 

n n 

shortest single-erasure-error-correcting code, the ((4,2^,2 )) c ode |9|, llO|, and the shortest 



single-arbitrary-error-correcting code, the ((5,2^,3)) code H 



I2I. 



However, the restriction that the code dimension K of stabilizer codes is always a power 
of two indicates that these codes might not be optimal in many cases. In 1997, a nonad- 
ditive code, i .e., a code outside the class of stabilizer codes, with parameters ((5,6,2)) was 
constructed 13|. For the given parameters n = 5 and d = 2, this code outperforms the 
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best possible stabilizer code whose dimension is K = 2'^ = 4, providing the first example 
that good nonadditive codes exist. In 2007, after ten year of respite, other good nonaddi- 
tive codes were found, including a family of single-erasure-error-correcting codes containing 
a code ((5,5,2)) [l^, and a ((9,12,3)) single-arbitrary-error-correcting code [l^. In late 
2007, a general framework for constructing nonadditive codes, namely, the codeword stabi- 
lized (CWS) codes, was developed 16|]. These codes encompass stabilizer codes as well as 
all good known nonadditive codes. In addition, the CWS framework provides a powerful 
method to construct good nonadditive QECCs in a systematical way, and many good codes 



outperforming the best possible stabilizer codes were found |l6l-l22|. 

When implementing general QECCs, one big challenge comes from the complexity of the 
quantum circuits, for both encoding and decoding. To verify and enhance the control ability 
to implement complex quantum circuits are hence critical tasks for implementing QECCs in 
building scalable quantum computers. Coherent control of the simplest single-erasure-error- 
correcting code, the ((4,2^,2)) stabilizer code, has been demonstrated in optical systems 
231]. Coherent control of the simplest single-arbitrary-error-correcting code, the ((5,2\3)) 
stabilizer code, has been implemented in nuclear magnetic resonance (NMR) systems 24 1. 

Here we implement a ((5,5,2)) CWS quantum code to benchmark the ability of coher- 
ence control in NMR systems. Experiments are performed in a seven-qubit system where 
five spins are used to represent the five qubits of the code. This ((5, 5, 2)) code is one of the 
simplest nonadditive codes which encode a larger dimension of the Hilbert space than the 
best possible stabilizer mentioned above, for the given parameters n = 5 and d = 2, the best 
possible stabilizer code has dimension K = 2"^ = 4. Our experiment thus gives the first ex- 
perimental demonstration of a nonadditive code that outperforms the best possible stabilizer 
codes. Implementing one round of experiment typically requires approximately 320-440 ra- 
dio frequency (r.f.) pulses with a total duration of approximately 0.55-0. 72s, which is much 
longer than for the ((5, 2^, 3)) experiment with 368 r.f. pulses amounting to a total duration 



of 0.38s 



241 ]. Given that the decoherence time T2 of our NMR system ranges from about 



0.80s to 1.1s, our experiment is quite challenging with respect to coherence control of the 
system. Surprisingly, although only about 10-20% signal strength remains available, there 
is still clear evidence of quantum coherence, which is robust against the artificial errors 
generated by single-qubit rotations. What is more, we can partially characterize the errors. 
Our results demonstrate the ability to control NMR spin systems for implementing complex 
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quantum circuits for quantum computing [25|. 



QUANTUM ERROR-CORRECTING CODES 



The ((5,5,2)) code we have implemented has been first introduced in [IJ], and shown 



Q. 



to be a CWS code in [16|. This code uses five qubits to protect an arbitrary state in a 
five-dimensional Hilbert space, and is capable of correcting a single erasure error. In other 
words, if any of the five qubits is lost, we can always recover the arbitrary state in the five- 
dimensional code space as long as we know which qubit was lost. This equivalently means 
that this code corrects an arbitrary error on any of the five qubits given the error position 
is known. 

The entire procedure to demonstrate has three parts: encoding the input state into the 
five-qubit code, an arbitrary error on one of the five qubits at a known position, and decoding 
to recover the input state. The outline of the procedure is given in Figure [H We assign 
qubits 2-4 as register qubits for carrying the five-dimensional quantum state to be protected, 
and choose the five basis states as {|&)}b=0' which in binary form are 



{|000), |001), |010), |011), |100)}. 



(1) 



Then the input state |0in)234 can be an arbitrary superposition of these five basis states. 
Qubits 1 and 5 are syndrome qubits, which are initially both in the state |0). 

The encoding procedure maps the basis states {|&)}b=o ^'^ ilVb)} 
of the code. Here 



6=0; which form a basis 



(see also 
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but note the re-ordering of the code basis and the different input states for 



the encoding circuit). Using the general method to obtain the encoding circuit for a CWS 
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code [16| together with some optimization, for the ((5,5,2)) code we get the circuit shown 
in Figure [21 

For each location of the error (on one of the five qubits), a different decoding circuit is 
needed. These five decoding circuits are shown in Figure [3l for correcting errors happening 
on qubits 1-5, from top to bottom, respectively. The error operations, denoted as are 
also shown in the decoding circuits, separated by the double vertical lines, noting that £ 
can be an arbitrary unitary as well as non-unitary operation. 

We consider arbitrary unitary errors which can be represented as 

£ = e^'Rnie). (3) 

Here Rn{9) = e"*^" *^/^ denoting a rotation by 9 along the h axis, where a denotes a vector 
with three components formed by the Pauli matrices X, Y and Z. 
When the input state is chosen as 

|^6) = |0)l|6)234|0)5, (4) 

the output state after decoding has the form 

1*6) = |<f)l5|fe)234, (5) 

which recovers the input state \b) of qubits 2-4. The syndrome state |$)i5 contains infor- 
mation on the error. 

Expanding the unitary £ given in Eq. ([3]) as a linear combination of the matrices E, X, 
Z, and Y gives 

£ = cqqE + coiX + ciqZ + ciiF. (6) 
Here E is the identity operator acting on a qubit, and 

Coo = e'" cos(e/2), coi = -ze^" sin(e/2)n^, 

Cio = — «e*" sin(6'/2)n2, and cn = — ie*" sin(6'/2)ny, 

where fi = {nx,nx,nz). We then obtain the syndrome part |$)i5 as 

I$)i5 = CoolOO) + CoilOl) + CiollO) + cnlll). (7) 
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EXPERIMENTAL PROTOCOL 



Our experiment uses a Bruker DRX 700 MHz spectrometer. We choose ^^C-labelled 
trans-crotonic acid dissolved in d6-acetone as the qubit system [26|, l27|. The structure of the 
molecule is shown as the inset in Figure HI The methyl group (denoted as M) can be treated 
as a spin half nucleus using a gradient-based subspace selection 26|. The Hamiltonian of 
the seven spins can be represented as {h = 1) 



i<j 



(8) 



where Ui denotes the chemical shift of spin i, and Jij denotes the scalar coupling strength 
between spins i and j. The values of z/j, Jjj, and relaxation times are listed in Figure |H We 
exploit the seven spins as seven qubits in experiment. 



We prepare a labelled pseudopure state ps = OOZOOOO using the method in Ref. |26 |. 
Here the labelled qubit is in state Z and the other qubits are in state = |0)(0|. From left 
to right, the order of the spins is M, Hi, H2, Ci, C2, C3, C4. One should note that we are 
using the deviation density matrix formalism 28[. In experiment we do not use spins Hi and 
H2 in implementing the QECC, after the preparation of the state p^- These two spins are 
only affected by the hard refocusing pulses. We therefore concentrate on the subsystem of 



the five spins M, C1-C4, denoted as qubits 



some basic results on quantum circuits 29 



-5, in which the QECC is implemented. Using 



30[, we transform the circuits for implementing 



the QECC shown in Figures [2] and [3] into pulse sequences, shown in Figures [9] and [lO] in the 
Appendix for encoding and decoding, respectively. The single-spin rotations along 2;-axis 
are implemented through the evolution of the chemical shifts in the Hamiltonian of the 



spin system I3I 
r.f. pulses [31 



. Rotations along x- and y-axes are implemented by standard Isech-shaped 
33| for spins M and Ci, and pulses generated by the gradient ascent pulse 



engineering (GRAPE) algorithm 34] for C2-C4. The evolutions of the J-couplings between 
the neighboring qubits are implemented by numerically optimized refocusing pulses, which 
are standard Isech and Hermite-shaped r.f. pulses [3l|, l35|] and are not shown here. We 



choose artificial errors 



24 



to demonstrate the QECC, i.e., the error operations £ in Eq. 
(|3]) are generated by r.f. pulses. All operations are combined in a custom-built software 
compiler to minimize error accumulation 



27 



We protect the quantum coherence to demonstrate the QECC by putting one qubit of 
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the input state in a superposition |+) = (|0) + \l))/\/2. The corresponding coherence can 
be observed directly in the spectra of the input state and the final state after the completion 
of QECC (see below). We avoid the full quantum state tomography, which would make 
the demonstration more difficult because of the required additional experiments and the 
limitation of the spectral resolution of certain spins. We choose three input states as 

lV'.) = |0)|sfc)|0) (9) 

with 

|si) = |+)|00), |S2) = |01)|+), |S3) = |00)|+), (10) 

noting that all the basis states in Eq. ([T]) are included here. We introduce error operations 
of Y- and Z-type by setting n in Eq. ([3]) along the x-, y-, and 2;-axes, respectively, and 
choose a = 0. Using Eq. ([5]) one can obtain the output states represented as 

l^fc) = [C0s(e/2)|00)i5 - 2Sin(^/2)|j7)i5] ® \Sk)2M. (11) 

where = |01), |10), and |11) for X-, Z-, and F-type errors, respectively. Obviously one 
can obtain the results for the case of no error [E-tjpe error) by setting 6' = in the above 
equation. 

In NMR experiments, we can detect the states of the qubits through measuring the bulk 



magnetization 



37|, 



M(t) oc Tr[p(t) ^{Xj + zFj)e-*/^2,,]. (12) 



Here p{t) = U{t)p{0)W (t) where U{t) = e"**-^^^^^, and p(0) denotes the density matrix 
of the state to be measured after the r.f. pulses are switched off. Tg*^- denotes the effective 
transverse relaxation time of spin j. One can obtain the NMR spectrum through the Fourier 
transform of M{t). From some calculations, one can find that the single coherence elements 
in p(0) can be directly measured in the spectrum, and the peaks of a certain spin are 
associated with the states of the other spins because of the J-couplings. 

We first take spectra of Ci for input \ipi), and spectra of C3 for \ip2) and \ipz)) respectively. 
There is only one peak in each spectrum, coming from the state |0)|+)|000), |001)|+)|0), or 
|000|+)|0), in the input state \ipi), li^i), or I'^s), respectively, illustrated by the spectrum 
shown as the solid curve in Figure [5] (a) obtained from \ip2)- The amplitude of the peak is 
proportional to the element of the quantum coherence to be protected in state \sk) in Eq. 
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( ITO|) against the error operations E. We take the amphtude of the peak in the spectrum 
obtained from the input state as the reference to normahze the corresponding output signals 
after completion of the QECC procedure. 

In principle, full quantum process tomography (QPT) is required to completely evaluate 
the implementation of the QECC and analyze experimental errors. Based on the result 
of the QPT, one can estimate the experimental errors of the quantum gates for imple- 
menting the QECC through proper models for relaxation effects [ij. In practice, however, 
this task is rather difficult because the number of experiments required for the QPT in- 
creases exponentially with the involved qubits. In our five qubit case, more than one million 



(4^" — 4" = 4^° — 4^) real numbers have to be determined 



Moreover, the limitation of 



the spectral resolution of certain spins in our molecule further increases the difficulty. More 



recent strategies for characterizing quantum processes [39|, |40| might allow an improved 
evaluation of our QEC procedure. Here we choose a simplified approach to demonstrate 
robustness of the scheme against artificially introduced errors, and estimate experimental 
errors in combination with simulation results. 

Noting that in Eq. ( ITT]) the two syndrome qubits are not entangled with the register 
qubits, one can extract the state of the register by tracing over the syndrome qubits in 
In experiment we only take the spectra of the register qubits, and implement the 
partial trace by adding the splitting of the peaks of the register qubits caused the coupling 
of the syndrome qubits. Through Eq. (Ill]), one finds that the state |0)|+)|000), |001)|+)|0), 
or |000)|+)|0) contributes a peak with amplitude 

Ao = cos2(e/2), (13) 

and |j)|+)|00/), |j)|01)|+)|/), or |jOO)|+)|/) contributes a peak with amplitude 

= sin2(^/2) (14) 

in the state I^E'i), 1^2); or l^s); respectively, after the completion of the QECC, where the 
sum y4o + v4i = 1 shows the robustness against the error operations. We can therefore 
benchmark the QECC procedure through measuring Aq and A\. Moreover one can obtain 
the rotation angle Q through and A\ up to its sign. One should notice that we do 
not measure the syndrome qubits directly. The NMR spectrum of the output state is not 
necessary exactly in phase with the reference spectrum, due to experimental errors, e.g. 
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imperfection of r.f. pulses. We therefore choose the absolute values Iq = \Aq\, Ii = \Ai\, 
and I = \Aq + All represent the results. 

EXPERIMENTAL RESULTS 

We summarize our experimental results in this section. We have implemented the 
((5,5,2)) code in three different situations, for different input states and artificial errors. 

Setting A. A single input state \ip2) as given in Eq. and Pauli errors X, Y, Z, and 
E on each of the five qubits. The results are summarized in subsection A. 

Setting B. A single input state \ip2) as given in Eq. ([9]), and arbitrary X-, Y-, and Z-type 
errors, i.e., arbitrary rotations by an angle 6 along the x-, y-, and 2;-axes. We analyze the 
information on the rotation angle 6 and sumamrize the results in subsection B. 

Setting C. Three input states \ipi), \ip2), I^As) as given in Eq. ([9]), and arbitrary rotation 
errors along the y-axis. The results are summarized in subsection C. 

For all three situations, our results clearly demonstrate quantum coherence for imple- 
menting the QECC procedure. 

A. Correcting Pauli errors 

Figure [5] shows the experimental NMR spectra of C3 to illustrate the results for correcting 
Pauli errors, i.e., E, Z, X, and Y errors. In Figured (a) the spectrum shown as the thin curve 
is the reference spectrum, obtained from the input state \ip2) = |0)|01)|+)|0). The spectrum 
shown as the thick curve is obtained from the state of the equal weight superposition of all 
the computational basis states (i.e. |+)|+)|+)|+)|+)), where a scale factor of 16 is applied 
for better visualization. One can exploit the positions of the peaks in this spectrum for 
facilitating to locate the peaks in other spectra. Figures [5] (b)-(e) illustrate the results 
of error correction for E, Z, X, and Y errors, respectively. The shifts of the main peaks 
characterize the type of errors. The distance of the shifts, measured in J-couplings, is 
indicated by the arrows. 
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Error location 


1 


2 


3 


4 


5 


ao (experiment) 


0.156 ±0.020 


0.160 ±0.038 


0.191 ±0.036 


0.177 ±0.010 


0.142 ±0.038 


Qi (experiment) 


0.168 ±0.023 


0.176 ±0.019 


0.215 ±0.025 


0.177 ±0.038 


0.099 ±0.016 


/ (experiment) 


0.146 ±0.012 


0.164 ±0.018 


0.193 ±0.028 


0.174 ±0.027 


0.107 ±0.029 


ao (simulation) 


0.288 ± 0.031 


0.312 ±0.018 


0.230 ±0.016 


0.322 ±0.019 


0.249 ±0.015 


ai (simulation) 


0.295 ± 0.037 


0.369 ± 0.022 


0.218 ±0.020 


0.295 ± 0.016 


0.250 ±0.013 


/ (simulation) 


0.304 ± 0.007 


0.337 ±0.018 


0.222 ± 0.003 


0.315 ±0.010 


0.246 ± 0.006 



TABLE I: The fitted oq, cti and / in experiment and by simulation in case B. 



B. Correcting X-, Y-, and Z-type errors 

The input state is chosen as \ip2)- The results for correcting X-, Y-, and Z-tjpe errors 
with arbitrary rotations are shown in Figure [6l We exploit the averaged Jq, h and /, 
respectively denoted as Jq, /i and /, over the results of correcting the three types of errors 
as benchmark. The experimental results for correcting errors on qubits 1-5 are shown as 
Figures [6] (a)-(e), respectively, where Jq, /i, and I are indicated by squares, diamonds and 
circles. In comparison with the theoretical values Jg^ and J*^, the experimental data can be 
fitted as ao-^o'^ and aill'^ shown as the dash-dotted and dashed curves. We fit the data for 
/ using a constant function, shown as the plane lines. The remaining quantum coherence is 
robust against the rotation angle 9, demonstrating the success of the QECC procedure. The 
values of the fitted a^, ai and / are listed in Tabled We also list the results by simulation 
in Figures [6] (f)-(j), for correcting the errors happening on qubits 1-5, respectively. In 
simulation, we take into account the T2 effects, imperfection in refocusing protocols and the 
theoretical infidelity of the numerically generated pulses. The fitting results for ao, ai and 
/ are also listed in Table [B 



C. Correcting y-type errors from various input states 

We measure Jq, /i, and / for the three input states with k = 1, 2, 3, and exploit the 
average over the input states, respectively denoted as Jq, Ji and /, to represent the results 
of the error correction. The results in the experimental implementation and by simulation 
are shown in Figures [7] (a)-(e), and (f)-(j), corresponding to errors happening on qubits 
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Error location 


1 


2 


3 


4 


5 


ao (experiment) 


0.175 ±0.013 


0.130 ±0.011 


0.136 ±0.012 


0.161 ±0.022 


0.137 ±0.026 


Qi (experiment) 


0.165 ±0.017 


0.140 ±0.010 


0.128 ±0.010 


0.188 ±0.063 


0.093 ±0.019 


/ (experiment) 


0.155 ±0.012 


0.134 ±0.007 


0.119 ±0.018 


0.165 ±0.027 


0.108 ±0.024 


ao (simulation) 


0.314 ± 0.008 


0.330 ± 0.002 


0.276 ± 0.005 


0.324 ± 0.008 


0.269 ±0.031 


ai (simulation) 


0.301 ±0.012 


0.327 ±0.005 


0.243 ± 0.005 


0.311 ±0.017 


0.260 ±0.016 


/ (simulation) 


0.307 ± 0.008 


0.323 ± 0.004 


0.257 ±0.014 


0.314 ±0.008 


0.258 ±0.011 



TABLE II: The fitted oq, cti and / in experiment and by simulation in case C. 



1-5, respectively, where /q, h and / are indicated by squares, diamonds and circles. In 
comparison with the theoretical values Iq^ and 7*'^, the experimental data can be fitted as 
OqIq^ and ailf^ shown as the dash-dotted and dashed curves in Figures [7] (a)-(e). The data 
for I can be fitted as a constant function, and the fitting results are shown as the plane lines. 
The remaining quantum coherence is robust against the rotation angle 6, demonstrating the 
success of the QECC procedure. The values of the fitted ao, ai and I for the results in 
experiment and by simulation are listed in Table Ull 

DISCUSSION OF THE EXPERIMENTAL RESULTS 

A. Extracting partial information of error operations 

Through measuring the qubits of the register [see Eqs. ( !T3|) and (fT4|) ]. we can obtain 
partial information of the errors, i.e., the types of the errors illustrated in Figure O and the 
absolute values of the rotation angles of the errors, exploiting the couplings between the 
register and syndrome qubits. Figure |8] illustrates the experimental results obtained from 
the data in Figure [7] for extracting absolute values of the rotation angles. In Figure [H] (a), 
the theoretical expectation = 6' is shown as the line. The experimental data is distributed 
near the line. We fit the data using a linear function as Q = a6 + b. The coefficients a 
and b are shown in Figure [8] (b) and (c), respectively, where the dashed line indicates the 
theoretical value of a = 1. The fitting results show a good agreement between experiment 
and theory. 
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B. Error sources in implementation 

The simulation is mainly used to analyze experimental errors. Excluding the preparation 
of the pseudopure state, the number of r.f. pulses ranges from about 320 to 440, and the 
duration of the experiments ranges from about 0.55s to 0.72s. The duration is comparable 
with, and indeed very close to the transverse relaxation times (T2's in Figure |4]). Hence the 
limitation of the coherence time contributes major experimental errors. Through comparing 
the results by simulation with and without T2 effects, we estimate that the limitation of 
T2 contributes to about 34%-45% of the loss of signal. Additionally, the errors due to 
imperfections of the refocusing protocols and the implementation of r.f. pulses contribute 
about 10%-16% and 3%-19% of the loss of signal, respectively. 



C. Comparison with the implementation of the ((5,2\3)) code 
We use the same seven- spin- qubit NMR system that was used in the implementation 



of the ((5,2^,3)) codes 2J]. It should be emphasized that in the ((5,2^,3)) experiment, 
the entire circuit one would need to carry out consists of four parts: encoding, (artificial) 
errors, decoding, and recovery. The decoding circuit is just the reverse of the encoding 
circuit, but the recovery circuit involves approximately three times more gates than the 
encoding/decoding circuit. However, because the artificial errors demonstrated for the code 
((5, 2^, 3)) are all Pauli errors, after the decoding circuit using only Clifford gates (gates that 
map Pauli operators to Pauli operators), the five qubits are in a product state. Only a single 
qubit carries information that needs a recovery operation, while the other four qubits carry 
error syndrome information. Consequently, the noise (mainly dephasing noise) during the 
recovery procedure does not have much effect on the final experimental fidelity, as coherence 
is only needed to be maintained on a single qubit. 

On the contrary, we use three qubits to carry information. Although our experiment 
requires implementing quantum circuits some of which are of comparable size as those for 
the ((5,2^,3)) code, quantum coherence in a larger Hilbert space has to be maintained 
throughout the entire procedure of the code. In implementation, maintaining coherence in 
a larger space e.g. three qubits, is more difficult than in a one-qubit space, because the high 
order coherence could decay faster than the single coherence under effects of decoherence |41 |. 
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This case might explain the low remaining signals in our experiment (10-20%), compared 
with the ((5, 2^, 3)) experiment where the intensity of the remaining signals ranges 48-87%. 

CONCLUSION 

We implement a five-qubit quantum error-correcting code using NMR. The code protects 
an arbitrary state in a five-dimensional Hilbert space, and is capable of correcting a single 
erasure error. The code is a CWS code which encodes a larger Hilbert space than any 
stabilizer code with the same length and being capable of correcting a single erasure error. 
Our results demonstrate a good nonadditive quantum code in experiment for the first time. 

Compared with the previous implementation of another five-qubit code using the same 



seven-spin-qubit system 2J], which encodes a two-dimensional Hilbert space and is capable 
of correcting an arbitrary single-qubit error, the pulse sequences in our experiment are 
more complex. In order to shorten the length of sequences, we exploit pulses optimized by 
the GRAPE algorithm to implement the 7r/2 spin-selective pulses for certain spins. The 
duration of our experiments ranges from 0.55s to 0.72s, which challenges the ultimate limit 
of coherent control of the system, with a typical T2 time ranging from 0.8s to 1.1s. Despite 
experimental imperfections which induce signal loss, the signal remains resilient against the 
artificial errors. 
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Appendix: Pulse sequences 

We list the pulse sequences for implementing the encoding and decoding circuits as Figures 
iQlandfTOl respectively. 
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|0)i 

I (/'in) 234 < 
|0)5 

FIG. 1: Outline of the protocol for the QECC procedure. Five qubits 1-5 are represented by five 
lines, from top to bottom, where qubits 2, 3, and 4 are exploited to carry the input state |</>in)234) 
and qubits 1 and 5 are syndrome qubits. The error can happen on any of the five qubits. The 
decoding operation depends on the location of the error. In the ideal case the output state obeys 

\^ova) = \(t>hi)- 



encoding 



error 



decoding 




FIG. 2: Quantum circuit for encoding. Elements © and • connected by a line denote a controlled- 
NOT or Toffoli gate conditioned on the state |1). X denotes the Pauli matrix (Jx-, acting as a NOT 
gate, and H denotes a Hadmard gate. 
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FIG. 3: From top to bottom, the quantum circuits on the right side of the double vertical lines are 
for decoding and correcting the errors happening on qubits 1-5, respectively. £ on the left side of 
the double vertical lines denotes the error operation. 
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FIG. 4: (Color online) Characteristics of the molecule of trans-crotonic acid. The molecular struc- 
ture is shown as the inset. The chemical shifts and J-coupling constants (in Hz) are on and above 
the diagonal in the table. The longitudinal and transversal relaxation times Ti and T2, which are 
listed at the bottom, are estimated by the standard inversion recovery and Hahn echo sequences. 
The chemical shifts are given with respect to reference frequencies of 700.13 MHz (protons) and 
176.05 MHz (carbons). The nine weakly coupled spin half nuclei can provide seven qubits since 



the methyl group can be treated as a single qubit using a gradient-based subspace selection 
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FIG. 5: (Color online) Experimental NMR spectra of C3 to illustrate the results for correcting 
Pauli errors. In figure (a) the spectrum shown as the thin curve is obtained from the input state 
|'02) = |0)|01)l+)|0)- The spectrum shown as the thick curve is obtained from the state of the equal 
weight superposition of all the computational basis states, with an applied scale factor of 16 for 
better visualization. Figures (b)-(e) illustrate the results of error correction for E, Z, X, and Y 
errors, respectively. The shifts of the main peaks, which are indicated by the arrows, characterize 
the type of errors. 
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FIG. 6: (Color online) Averaged results of the correction for X-, Y- and Z-type errors. The input 
state is chosen as IV'2)) and the signal in the input state is chosen as reference to normalize the 
signals after the completion of the QECC procedure. The left and right columns show the results 
in experiment and by simulation, respectively. The results for errors on qubits 1-5 are shown from 
top to bottom, respectively. The averaged Iq and /i, denoted as Iq and Ji, are indicated by squares 
and diamonds, respectively. In comparison with the theoretical values /g'* and 7*^, the data can 
be fitted as oqIq^ and ai/*^ shown as the dash-dotted and dashed curves in figures (a)-(e) for 
experiment and (f)-(j) for simulation, where ao and ai are constant coefficients in each figure and 
listed in the main text. The averaged /, denoted as / and marked as circles, indicates the intensity 
of the remaining signals after completion of the QECC procedure. The remaining signal is robust 
against the rotation angle 9, demonstrating the success of the QEC procedure. 
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FIG. 7: (Color online) Averaged results over various input states {ipk) {k = 1,2,3) for correcting 
y-type errors. The left and right columns show the results in experiment and by simulation, 
respectively. The results for errors on qubits 1-5 are shown from top to bottom, respectively. We 
choose the signals of the input states as reference to normalize the signals after the completion of the 
QECC procedure. Iq and Ii are indicated by squares and diamonds, respectively. In comparison 
with the theoretical values Iq^ and if^, the data can be fitted as oqIq^ and ailf^ shown as the 
dash-dotted and dashed curves in figures (a)-(e) for experiment and (f)-(j) for simulation, where 
ao a.nd ai are listed in the main text. The intensity / of the remaining signal, marked as circles, 
is robust against the rotation angle 6, demonstrating the success of the QEC procedure. 
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FIG. 8: (Color online) Experimental results of error syndromes. Using the measured Jq and Ii, 
we obtain the absolute values of the rotation angles for the errors, denoted by 0. The theoretical 
expectation Q = 6 is shown as the line. The experimental data for errors happening on qubits 
1-5 are denoted as filled circles, triangles, diamonds, stars, and empty circles, respectively. The 
deviation of the data for errors happening at qubit 5 is larger than the other cases, because more 
r.f. pulses are required for correcting the errors at qubit 5 (see Figure [T0|) . We exploit Q = a9 + b 
to fit the data. The coefficients a and b are shown in figure (b) and (c), respectively, where the 
dashed line indicates the theoretical value o = 1. 
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b) Error at qubit 2. 
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c) Error at qubit 3. 
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d) Error at qubit 4. 
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e) Error at qubit 5. 
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FIG. 10: Pulse sequences to decode the ((5,5,2)) code, for errors happening on each of the five 
qubits. 



